Ultra-fast spin dynamics: the effect of colored noise. 
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Recent experimental results have pushed the limits of magnetization dynamics to pico- and fem- 
tosecond timescales. This ultra-fast spin dynamics occurs in extreme conditions of strong and rapidly 
varying fields and high temperatures. This situation requires new description of magnetization dy- 
namics, even on a phenomenological level of the atomistic Landau-Lifshitz-Gilbert equation, taking 
into account that the correlation time for electron system could be of the order of the inverse charac- 
teristic spin frequency. For this case we introduce the thermodynamically correct phenomenological 
approach for spin dynamics based on the Landau-Lifshitz-Miyasaki-Seki equation. The influence of 
the noise correlation time on longitudinal and transverse magnetization relaxation is investigated. 
We also demonstrate the effect of the noise correlation time on demagnetisation rate of different 
materials during laser-induced dynamics. 

PACS numbers: 75.40. Gb; 75.40.Mg 



One of the fundamental questions of modern solid state 
physics is how rapidly the magnetization can respond to 
an external excitation. The recent development of time- 
resolved pump-probe experimental techniques using X- 
ray spectroscopy based on synchrotron radiation [l| and 
the Stanford linear accelerator (SLAC) 0,0] has allowed 
the investigation of magnetization dynamics on the pico- 
second time scale. The use of the powerful femto-second 
lasers 0, [|| 0| has pushed this limit down to the femto- 
second time-scale. The physical processes underlying the 
response of the magnetization on this ultra-short time- 
scale are complicated and far from being understood, 
but clearly involve the excitation and consequent non- 
equilibrium interaction of electron, phonon and spin sub- 
systems. Spin dynamic processes on this time-scale oc- 
cur under extreme conditions remarkably different from 
those typical for dynamics at longer time-scales. Firstly, 
the three subsystems (electron, phonon and spin) are 
not in equilibrium with each other. Secondly, spin dy- 
namic processes occur under very strong fields with dif- 
ferent sources. In particular, in the experiments using 
the SLAC the magnitude of the external field can be as 
large as 20T. The strongly non-homogeneous magnetiza- 
tion processes are driven by the exchange field, having a 
magnitude > 100T. Finally, in the laser-induced magne- 
tization dynamics the effective temperature is increased 
up to and often above the Curie temperature [1, 0| . 

Atomistic spin models have proved to be a power- 
ful approach to model ultra-fast magnetization dynamics 
For example, in the case of laser-induced magneti- 
zation changes, spin models provide important physical 
insight into the spin- reordering process, establishing the 
linear character of the demagnetization during the sub- 
picosecond regime and predicting the origin of different 
recovery rates in the pico-second regime. The basis of 
these models is the stochastic Landau-Lifshitz (LL) equa- 



tion for each localized magnetic moment s^: 
Si = 7[s f x - ja [sj x [s, x H 2 



(1) 



Here is the local effective field which includes Zee- 
man, exchange, anisotropy and magnetostatic contribu- 
tions, augmented with a stochastic field with the 
following properties for both components and different 
spin sites: 

&(*)>= 0, <&(t)&(0> = ^^(<-*')*«- (2) 

Here T is the temperature, 7 is the gyro-magnetic ra- 
tio, jjL s is the magnetic moment and a is the parameter 
describing the coupling to the bath system. The basis 
of this equation is the separation of timescales, assum- 
ing that the bath (phonon or electron system) is much 
faster than the spin system. In this case, the bath de- 
grees of freedom can be averaged out and replaced by 
a stochastic field with white noise correlation functions. 
The coefficient in front of the delta function in Eq.(f2]) is 
determined by the fluctuation-dissipation theorem. The 
assumption of white noise is therefore invalid for magne- 
tization dynamics occurring on a timescale comparable 
to the the relaxation time of the electron system. The 
typical correlation time for the electron system in metals 
is ^ 10/ s [Hi- Such magnetization dynamics time scale is 
now commonly achieved by applying femtosecond laser 
pulses. A further limitation of this approach comes from 
the fact that characteristic frequencies of the magnetiza- 
tion process are now also of the order of the timescale 
of the noise (electron) variable. Therefore, for modelling 
the ultra-fast magnetization experiments the approach 
(J]) could break down. 

The aim of the present Letter is to present a classical 
formalism beyond the white noise approximation to form 
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a strong physical basis for models of ultra-fast magne- 
tization dynamics in extreme conditions. First, we in- 
troduce the formalism and show that our approach is 
consistent with the equilibrium Boltzmann distribution 
and coincides with the previous atomistic approach for 
small correlation times for the bath variable. As main 
implications of the correlated noise approach, we discuss 
the influence of noise correlations on the most relevant 
characteristics of magnetization dynamics: the longitu- 
dinal and transverse relaxation times. Finally, we model 
the laser-induced demagnetisation rate for materials with 
different noise correlation times. 

The standard generalization of the white noise to in- 
clude correlations is the Ornstein-Uhlenbeck stochastic 
process [9]. However, when we implemented this pro- 
cess within the LL dynamics, we have seen that in agree- 
ment with the general theory [IcJ, the Boltzmann distri- 
bution at equilibrium is not recovered for magnetization 
dynamics with correlated noise in this case. The devi- 
ations invariably correspond to precessional frequencies 
of the order of the inverse correlation time. The colored 
noise approach based on the Ornstein-Uhlenbeck process 
may provide a reasonable description in some situations 
when the system goes to a stationary condition not nec- 
essary coinciding with the equilibrium one (an example 
of this could be the spin-torque pumping into a magnetic 
system). However, in experiments such as those corre- 
sponding to laser-pulsed induced dynamics, a stochastic 
approach giving the correct equilibrium magnetization 
long after the laser pulse is gone, is necessary. 

A suitable approach has been found in the work of 
Miyazaki and Seki [111 who generalized the Langevin 



equation for one spin to a non-Markovian case. The ap- 
proach has been introduced for one spin at high tem- 
peratures, neglecting the interactions with other spins 
and assuming that their role is to provide the bath en- 
vironment. In the present paper we generalize this ap- 
proach to a many spin case, similar to the standard way 
of Eq.JTJ where the applied field is substituted by the 
local field. We assume that the bath variable is due to 
external sources such as electrons. The other assumption 
made in this approach is that the spin is connected lo- 
cally to the bath. Consequently, the set of equations for 
magnetization dynamics (in the following called Landau- 
Lifshitz-Miyazaki-Seki (LLMS)) reads: 



= 7[sj x (Hi + r/J], 
Vi = lVi-XSi)+Ri 



(3) 



with the fluctuation-dissipation theorem for the 
bath variable: (Ri(*)) = ; (Ri(t)R 3 -(i')) = 
(2xksT /T c )8ij8{t — t'). The parameter x describes 
the coupling of the bath variable to the spin. The 
precession term in the first equation of the set ([3]) has 
the same form as in the Eq.Q. However, the damping is 



now described by the second equation in this set where 
also the bath variable adjusts to the direction of the 
spin due to the interaction with it. In the limit r c — * 
the stochastic LL equation (TJ is recovered This 
also provides a relation between the damping and the 
coupling constants as a = JXTc, giving a more precise 
physical sense to the LL damping constant at atomistic 
level. 

For integration of Eqs.((3]) the Heun integration scheme 
was modified for this special case. First of all we investi- 
gated the equilibrium properties for an ensemble of non- 
interacting spins. In all cases of large fields, temperatures 
and correlation times, the correct Boltzmann distribution 
is obtained at equilibrium (see inset in FigQ]). 
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FIG. 1: Equilibrium magnetization as a function of tempera- 
ture for a system of N = 32 3 interacting spins, integrating the 
LLMS equation with different correlation times and integrat- 
ing the LL equation. The inset shows distribution functions 
for non-interacting spin system modeled within the LLMS ap- 
proach for different values of the reduced field £ = ^sH/ksT 
{^Htc = 1.76) and correlation time r c = lOfs. The solid line 
in the inset represents the Boltzmann distribution. 

Now we turn to the multi-spin system. First of all, 
we prove that the stochastic Eqs.([3]) for a multi-spin sys- 
tem are consistent with the standard equilibrium prop- 
erties. For this purpose we use the formalism of the On- 
sager kinetic coefficient method applied in Ref.[H] for 
the LL multispin system U). The system |(3|) is linearized 
near equilibrium and represented in a general form of the 
Langevin equation: 



(It 

(*(*)) 



(4) 



0; (ri(0)rj(t)) = miS(f) 



Here stands for small deviations of the stochastic vari- 
ables Sj or r\i from their equilibrium values, X, repre- 
sent their thermodynamically conjugate variables and 
My = lij + Iji- F° r the spin variable we have: Xj = 
— (n s /kBT)Hj , where Hj is the internal field correspond- 
ing to a particular lattice site and spin component. Un- 
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like Eq.JT]), the first equation in Eqs.((3]) contains only a 
precessional term and, therefore, the corresponding ki- 
netic coefficients are antisymmetric in spin components, 
giving for this equation /Xy = 0. Taking into account the 
generalization of the internal energy to include the bath 
variable as ^({s,}, {%}) = F ({s t }) + J2M li^x) -V&i], 
where Fo({sj}) is the internal energy without the bath 
variable, the conjugate variable to the bath one is Xj = 
— x s j)/(kBTx)- Therefore, the corresponding ma- 
trix of the kinetic coefficients is diagonal and for the sec- 
ond equation we obtain fXij = (2kBTx/T c )Sij . Conse- 
quently, we have proven that, under the supposition of 
local coupling of the spin to the bath variables, the set 
of multi-spin equations is consistent with the equilib- 
rium properties. 

In our simulations for the multispin system we use a 
Heisenberg Hamiltonian on a cubic system of 32 3 mag- 
netic moments with nearest-neighbor interactions only, 
with fig = 1.45/is and the Curie temperature T c = 700K 
(ksTc w 1.44J). The coupling parameter x was chosen 
to give the LL damping parameter a = 0.01. In FigQ] 
we present calculations of the equilibrium magnetization 
as a function of temperature for spin systems with differ- 
ent values of the correlation times. Independence of the 
equilibrium properties on the correlation time, and the 
agreement with calculations using the LL equation with 
uncorrelated noise demonstrates our generalization of the 
LLMS equation to multi-spin systems. Consequently, the 
LLMS equation provides a basis for the phenomenolog- 
ical description of magnetization dynamics in extreme 
situations of high temperatures, and large and rapidly 
varying external fields. The advantage of the approach 
is also that the fluctuation-dissipation theorem is not ap- 
plied directly to the spin variable. Therefore, the bath 
variable (for example electrons) and the spin system need 
not be in equilibrium with each other. 

Next we discuss the most important implications of 
this new approach to the ultra- fast dynamics. It is known 
that during the excitation with spatially inhomogeneous 
fields in the Terahertz range 0, [aland also during laser- 
induced magnetization dynamics^, 0] , strong local disor- 
dering of the spin system occurs. The dynamics in these 
cases is governed by field or temperature excited high- 
frequency spin waves which are responsible for the effec- 
tive damping. The important dynamical feature then is 
the rate of magnetization recovery. During these pro- 
cesses two types of relaxation could be distinguished. 
The first one known as longitudinal relaxation is respon- 
sible for linear magnetization recovery, i.e. the magne- 
tization magnitude. During the laser-induced demag- 
netization, the longitudinal relaxation is responsible for 
the femtosecond demagnetization. The second one is the 
transverse relaxation when the magnetization vector re- 
laxes to the direction parallel to the internal field via 
magnetization precession. The longitudinal relaxation 
time increases with the temperature while the transverse 
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FIG. 2: Longitudinal relaxation time (normalized to the un- 
correlated case) as a function of the correlation time for vari- 
ous temperatures calculated within LLMS approach. The in- 
set shows longitudinal relaxation for various correlation times 
and T = 600K. 



relaxation time has minimum at T c [13J. 

To simulate the longitudinal relaxation, we start with 
the initial condition {sf } = 1.0 and observe the system 
to relax at given temperature T. The obtained relaxation 
curves are then fitted to exponential decay to extract the 
longitudinal relaxation rate. The longitudinal relaxation 
time, normalized to the uncorrelated case, is presented 
in Fig|2] as a function of the noise correlation time. The 
longitudinal time calculated by means of the LL approach 
([1]) is of the order of lOfs (rjf = 28fs at T = 300K). 
For correlation time r c < lfs the uncorrelated approach 
gives the same results as the LLMS one. However, one 
can see that r c ~ 10 — lOOfs gives a dramatic increment 
of the longitudinal relaxation time. The effect is less 
pronounced at higher temperature since in this case the 
temperature contributes to the loss of correlations. 

Next, we investigate the transverse relaxation in Fig. 
[3l The transverse relaxation time is defined by the mag- 
netisation precession and normally is much slower than 
the longitudinal one. For one spin t± =Tj_[1+ (loht c ) 2 ] 
1 11 ] - where luh is the field-dependent precessional fre- 
quency. Consequently, the influence of the correlation 



time on the transverse relaxation may be expected only 
for strong applied field for which ujh ~ r^ 1 and, thus, 
could be relevant for SLAC experiments. To show the 
influence of the noise correlation on the transverse relax- 
ation, we model the precessional dynamics at strong ap- 
plied field H = 24.85T (0.05 J/ /J,,). In this case, the spin 
system was first equilibrated at given temperature and 
applied field. After that the whole system was rotated 
to an angle 30° and the relaxation to the direction paral- 
lel to the applied field was observed. For this particular 
strong applied field, the correlation times t c > 100/s are 
necessary in order to see their influence on the transverse 
relaxation. 

To demonstrate the immediate physical consequence 
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FIG. 3: Transverse relaxation time (normalized to the uncor- 
related case) as a function of the correlation time for various 
temperatures calculated within LLMS approach. The insert 
shows transverse relaxation for two temperatures T = 5K and 
T = 600K and r c = lOfs. 
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FIG. 4: Laser-induced demagnetisation as a function of time 
for various noise correlation times r c modelled within LLMS 
and LL approaches. 



/-electrons in rare-earths which have strongly relativistic 
nature. 

In conclusion, the standard phenomenological ap- 
proach to model spin dynamics has been generalized to 
the non-Markovian case. This approach is necessary in 
the extreme situations of large characteristic magnon fre- 
quencies occurring during ultra-fast magnetization pro- 
cesses. The advantages of the new approach are the fol- 
lowing: (i) the memory (correlation) effects arising from 
the fact that the bath variable responds to the spin direc- 
tion are taken into account. This corresponds to the sit- 
uation when the bath variable is not in equilibrium with 
the spin system, (ii) the fluctuation-dissipation theorem 
is not applied to the spin systems as in the standard LL 
approach (iii) In equilibrium the Boltzmann distribution 
is recovered. The price for this new approach is the use 
of two phenomenological constants: the phenomenologi- 
cal damping parameter a for the LL approach is substi- 
tuted by two phenomenological parameters in the LLMS 
approach: the correlation time t c and the coupling con- 
stant X- Several processes may be important in deter- 
mining these constants, as for example, the spin-orbit 
coupling, momentum relaxation, scattering rate and de- 
phasing time of conduction electrons. As in the LL ap- 
proach, these parameters will be material-specific and 
their physical origins should be clarified on the basis 
of first-principle approaches. We have shown that the 
ultra-fast magnetization dynamics is strongly influenced 
by these parameters which stresses the necessity of first- 
principle models, capable to clarify their physical origins. 

The authors acknowledge the financial support from 
Seagate Technology, USA and from European COST-P19 
Action. 



of the correlated noise for a more concrete experimen- 
tal situation, we present in FigQ] the modelling re- 
sults for the laser-induced demagnetisation for various 
noise correlation times and constant damping parame- 
ter a = 0.05. We suppose that the magnetization dy- 
namics is produced by a Gaussian laser pulse with 50 
fs duration and the fluency 33. 7m J /cm 2 . Similar to 
Ref.Q, we assume that the photon energy is transfered 
to the electrons and lattice but the magnetization is di- 
rectly coupled to the electron temperature T e . The lat- 
ter is calculated within the two-temperature model [? 

with the electron and lattice specific heat constants 
C e = 700J/m 3 K 2 x T™ ax , Q = 3 x 10 6 J/m 3 K and 
the coupling constant G e i = 8 x 10 17 W/m 3 K. Our re- 
sults clearly show strong impact of the noise correlation 
time on the degree of demagnetisation during the laser- 
induced process. Namely, the materials with small r c are 
demagnetised easier. This could be true, for example, 
for d-electrons in metals with large scattering rate or for 
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